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The alignment, ordering and rotation of elongated granular particles was studied in shear flow. 
The time evolution of the orientation of a large number of particles was monitored in laboratory 
experiments by particle tracking using optical imaging and x-ray computed tomography. The exper- 
iments were complemented by discrete element simulations. The particles develop an orientational 
order. In the steady state the time and ensemble averaged direction of the main axis of the par- 
ticles encloses a small angle with the streamlines. This shear alignment angle is independent of 
the applied shear rate, and it decreases with increasing grain aspect ratio. At the grain level the 
steady state is characterized by a net rotation of the particles, as dictated by the shear flow. The 
distribution of particle rotational velocities was measured both in the steady state and also during 
the initial transients. The average rotation speed of particles with their long axis perpendicular to 
the shear alignment angle is larger, while shear aligned particles rotate slower. The ratio of this 
fast/slow rotation increases with particle aspect ratio. During the initial transient starting from an 
unaligned initial condition, particles having an orientation just beyond the shear alignment angle 
rotate opposite to the direction dictated by the shear flow. 



I. INTRODUCTION 

Granular flows involving elongated particles are very 
important in agriculture, industry or natural processes. 
In a shear flow the interparticle contact forces often lead 
to orientational ordering and alignment of the particles. 
While the basic properties of granular flows have been 
widely studied for rapid (collisional), quasistatic and in- 
termediate regimes [TM24). the flow properties of sam- 
ples consisting of elongated objects is not well under- 
stood. After early experiments on grass seeds [6] and 
rice grains [7], recent numerical studies focused on shear 
flows of ellipsoidal particles [H1IIH1 HO]- These latter stud- 
ies revealed the topology of the contact forces [TH1 for 
sheared ellipsoidal particles and also described how the 
stresses and the corresponding particle orientations de- 
pend on the aspect ratio of the ellipsoids [5] ■ In a recent 
work [35] we have demonstrated the basic properties of 
the shear induced order focusing on the time and ensem- 
ble averaged properties only. We compared these to var- 
ious systems including colloids, nanoparticles or nematic 
liquid crystals. We have also shown that the shear align- 
ment results in a considerable reduction of the effective 
friction of the material. 

In the present paper we explore the behavior of indi- 
vidual particles: we detect their rotation induced by the 
shear flow, and measure the average rotation speed as a 
function of particle orientation during stationary shear 
and also during the transient when shear is applied to 
an initially unoriented or misoriented system. We also 
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present the basic properties of the shear aligned state for 
an extended set of materials. Details of our experimental 
and numerical methods are given. 



II. EXPERIMENTAL SETUP 
A. Experimental geometry and materials 

The granular material was placed in a so called cylin- 
drical split bottom cell (see Fig. [I]). In this geometry a 
circular plate at the bottom of the container is rotated 
with a constant speed. If the granular layer is thin enough 
(compared to the size of the rotating plate) the middle 
part of the granular sample will rotate with the plate like 
a rigid body. The outer part of the sample (near the walls 
of the cylindrical container) remains stationary. The re- 
gion in between these undisturbed parts is continuously 
sheared and is called the shear zone (or shear band) . The 
shear zone is indicated with dark gray in Fig.Qa). Other 
properties of this configuration, like the precise geometry 
of the shear zone and its dependence on the filling height 
were described in previous studies [5T1 [35] HI] • 

In our experiments, various prolate particles were in- 
vestigated as shown in Fig. [ljd) . The length to diameter 
(or aspect-) ratios of the particles were L/d = 2.0, 3.4 
and 4.5 (rice grains), 1.4 and 3.5 (glass cylinders), and 
2.0, 3.3, 5.0 (wooden pegs). Two experimental methods 
were used to visualize and characterize the shear-induced 
alignment: (i) complete 3D imaging of the sample with x- 
ray computed tomography (CT) and (ii) optical imaging 
of the particles at the upper surface of the system with 
a digital camera. Pegs were investigated by x-ray CT, 
glass cylinders and rice particles by optical imaging. The 
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FIG. 1: (color online), (a) Illustration of the experimental 
setup. The granular material consisting of elongated grains 
is sheared as the circular bottom plate continuously rotates. 
Shear strain is localized in the region marked with dark gray, 
(b) Example image of the surface of the shear zone for glass 
cylinders with d — 1.9 mm and L/d — 3.5. The typical orien- 
tation of the particles encloses a small angle 9 with respect to 
the streamlines, i.e. the particles are slightly inclined towards 
the velocity gradient, (c) Shear rate 7 as a function of the 
radius r; also glass cylinders, L/d = 3.5. Bins with different 
average shear rate are marked, (d) Photographs of the parti- 
cles studied, with L/d — 2.0, 3.4, 4.5 (rice), as well as 1.4 and 
3.5 (glass cylinders), and 2.0, 3.3, 5.0 (pegs). 



diameter of the rotating plate was D p ps 80 ef, while the 
container diameter was D c ~ 140 d in both experiments. 
We used a typical filling height of H ps lid sa 0.14 D p . 



B. Particle detection by x-ray computed 
tomography 

The complete 3D arrangement of all particles in the 
shear zone was determined experimentally for three sam- 
ples (pegs with L/d = 2.0, 3.3 and 5.0) using x-ray CT 
|30j . For these measurements we used a medical angiog- 



raphy system, i.e., a rotational c-arm based x-ray device 
equipped with a fiat-panel detector [31] . The spatial res- 
olution of the device was between 1.5 pixel/mm and 2 
pixel/mm. The primary purpose of these experiments 
was to show that the order and alignment data extracted 
from surface particles are representative for the global 
order in the shear zone, with only slight variations. The 
particles were detected by applying a watershed algo- 
rithm to the recorded images. The particles' orientation 
was defined as the direction of the largest eigenvector 
of the individual moment of inertia tensors. The shear- 
induced orientational order is monitored by diagonalizing 
the symmetric traceless order tensor T: 
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where 3- n ^ is a unit vector along the long axis of parti- 
cle n, and the sum is over all N detected particles. The 
largest eigenvalue of T is the primary order parameter 
S. The shear-induced alignment is characterized by the 
average alignment angle 9 av , which is the angle of the 
corresponding principal axis of T measured with respect 
to the flow direction [see Fig. [jjb)]. Positive 9 corre- 
sponds to inclination towards the velocity gradient. The 
components of T were determined both near the surface 
(particles with centers in the top ss 1 cm layer, approx. 
2 particle diameters) and in the bulk (2 cm thick layer 
below), corresponding to about 10,000 and 20,000 parti- 
cle positions, respectively. The experiment is described 
in detail in I3TH. 



C. Optical particle detection 

For rice grains and glass cylinders, the top surface of 
the whole sample was monitored with a digital camera 
using a frame rate of 18 fps. As an illustration, a sec- 
tion of an image is shown in Fig. [ijb). The shear rate 
7 is not constant across the shear zone, but it changes 
with the radius r measured from the center of the cell. 
It is expressed as 7 = dv(r)/dr — v(r)/r, where v(r) 
is the tangential velocity of the material. The tangen- 
tial velocity v(r) was determined using a self developed 
PIV (Particle Image Velocimetry) code and the time av- 
eraged shear rate was calculated. This shear rate profile 
is presented in Fig. [ljc). The shear zone was divided 
into about 20 bands (widths 0.19 cm sa d) as illustrated 
in Fig. [l|c). The average shear rate was determined for 
each band. The orientation and in-plane length of colored 
tracer grains (6% of the sample) were detected, yielding 
about 60,000 particle positions in each band for each type 
of grain. This way, a single experiment with constant ro- 
tation speed n provides information about the shear rate 
dependence of the alignment |23j . In order to cover a 
wide range of shear rates, the experiments were repeated 
at three different values of f2. 
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D. Torque measurement 

In order to monitor changes in the mechanical resis- 
tance of the sample against shearing, we continuously 
measured the torque needed to maintain a constant ro- 
tation rate of the bottom plate during the evolution of 
the alignment. This has been done the following way: 
The driving motor was mounted on a rotating table. To 
keep the body of the motor fixed and let the plate rotate 
in the sample container, we use a load cell attached to 
the body of the motor. We measure the torque needed 
to hold the motor while the granular material is sheared 
with a constant speed by the rotating bottom plate. 



III. NUMERICAL SETUP 

Discrete Element Method (DEM) simulations were 
performed to complement the experiments. A modified 
version |24j of the general purpose molecular dynamics 
code LAMMPS [27] was used to generate stiff elongated 
particles by gluing spheres together. Long particles gen- 
erated such way have parts with concave curvature and 
thus are prone to interlock if their structure is identical. 
Therefore we prepared polydisperse samples by allowing 
the variation of the following quantities: the radius of the 
spheres by 10%, the length of the particles by 10%, and 
the overlap of the spheres between 30-70% of the radius. 
We also varied the number of spheres required to build 
up a long particle. 




FIG. 2: (color online). Example cut of the numerical setup. 
The wall consisting of larger spheres imposes constant pres- 
sure and stationary shear to the sample. The snapshot was 
taken in the steady state, where most of the particles are shear 
aligned. 

The dynamics of 700-1200 elongated particles was 
modeled in simple Couette flow, where the parallel walls 
were moved in the opposite directions as indicated in Fig. 
[2] by the arrows. The walls consisted of twice as large 



spheres as the average particle diameters. In the other 
two directions periodic boundary conditions were chosen. 
The particles were confined with constant pressure during 
shear. The generated shear rate was relatively uniform 
across the sample [2"5] , 



IV. RESULTS 
A. Stationary state 

As described in Section|nJ a sufficiently large set of par- 
ticle positions and corresponding orientations were de- 
tected using either optical imaging at the surface or x-ray 
CT in the whole sample. The distribution of the align- 
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FIG. 3: (color online). Distributions of the orientational an- 
gle 9 of the particles with respect to the streamlines for (a)-(c) 
cylinders with L/d = 5.0, 3.5, and 2.0, and (d)-(f) rice parti- 
cles with L/d — 4.5, 3.4, and 2.0, respectively. For panels (a) 
and (c) data obtained by x-ray CT, other data obtained by 
optical particle detection. The vertical blue line corresponds 
to the stationary average orientation, O av - 

ment angles 9 (measured in the shear plane with respect 
to the local direction of the streamlines) is presented for 
cylinders with length to diameter ratios of L/d = 5.0, 
3.5, and 2.0 and and rice with L/d = 4.5, 3.4, and 2.0 
in Figs.[3]ja)-(c) and Figs[3jd)-(f), respectively. The first 
thing to note is the aspect ratio dependence of the distri- 
butions. The graphs show that increasing L/d leads to 
a narrower distribution, i.e. a higher degree of shear in- 
duced ordering. The time and ensemble averaged particle 
orientation av was determined by fitting the data by a 
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Gaussian (with the constant offset subtracted), and is il- 
lustrated with a vertical (blue/gray) line in Fig. [3^a)-(f). 
Thus, asymptotically the average orientation encloses a 
small angle av with the streamlines, as it is illustrated in 
Fig. [fjb) for L/d = 3.5. The second interesting observa- 
tion is that for cylinders (L/d = 3.5 and 2.0) the distribu- 
tion is broader than for ellipses with similar aspect ratio. 
This can be understood by analyzing the shortest sam- 
ple (L/d = 1.4) in Fig. [4] Here, the distribution has two 
peaks, see the zoomed inset of Fig. |4|a). This suggests 
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with similar aspect ratio (compare Figs. |3^c) and|3^f)). 

Figure [3] also illustrates the typical scatter of the angu- 
lar distributions obtained by the two methods (optical or 
x-ray measurements) . Data for all rice samples as well as 
for cylinders with L/d = 3.5 were obtained by optical de- 
tection and the histograms presented were created using 
about 200,000 data points, while the data for pegs with 
L/d = 5.0 and 2.0 were collected from the surface region 
of the whole shear zone by x-ray CT (about 10,000 data 
points). Due to the smaller number of data points, the 
scatter of the distributions is larger in these latter cases. 

As explained before, the center O av of the orientational 
distribution was determined by fitting a Gaussian, and 
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FIG. 4: (color online), (a) Distribution of the orientational 
angle 8 of the particles with respect to the streamlines for 
cylinders with L/d — 1.4. Data obtained by optical particle 
detection. The inset in panel (d) is the same data on a differ- 
ent scale. The vertical blue line corresponds to the average 
orientation, Q av . (b) Illustration of typical (average) parti- 
cle alignment for longer (L/d = 3.5) and shorter L/d — 1.4 
cylinders. 

that these grains align with their longest extension (di- 
agonal) in the preferential direction. Since we determine 
the particle orientation by measuring the orientation of 
the symmetry axis, the distribution of the detected an- 
gles has two peaks, approx. 2 arctan(l/1.4) = 71° apart, 
which correspond to the two particle configurations illus- 
trated in Fig. |4|b) for L/d = 1.4. The effective average 
alignment angle in this case is between the angles cor- 
responding to the two peaks, i.e. around 6 av = 29°. 
Getting back to the case of longer cylinders, the above 
described effect could be one of the factors leading to a 
widening of the distributions (especially for the case of 
cylinders with sharp edges) compared to the case of rice 



FIG. 5: (color online), (a) The average orientational angle 
Oav of the particles with respect to the streamlines and (b) 
the order parameter S as a function of the shear rate 7. Data 
obtained by optical measurements for glass cylinders with 
length to width ratios L/d = 1.4 and 3.5. Both, av and 
S are independent of the shear rate. 

the order parameter S, defined as the largest eigenvalue 
of T, was calculated. Note that S is related to the width 
of the angular distribution a. As seen in Fig. [5] av and 
S (as well as tr, not displayed) are practically shear rate 
independent across more than three decades of 7. The 
slight decrease of S observed for two samples of rice [25] is 
not observed here for cylinders. The data points were col- 
lected using three different rotation rates (fi = 0.26, 0.52 
and 4.2 rad/s) of the bottom plate. The spatial variation 
of 7 in the shear zone allows to span a relatively wide 
range of 7 for each rotation rate, as marked in Fig. p[b). 
In terms of inertial number (defined as I — jdy p/P, 
where p is the density of particles and for the hydrostatic 
pressure P we use an indicative value taken at depth d), 
our measurements span the range between 0.0002 to 0.4. 
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The fact that the average orientation of the particles 
in the aligned state is shear rate independent is simi- 
lar to the well studied flow alignment of nematic liquid 
crystals. We have quantitatively compared our observa- 
tions to the properties of nematics in |25j . here we only 
wish to stress that the strong similarities with nemat- 
ics indicate that the shear alignment is robust and is of 
geometrical origin. Namely, even if the aspect ratios of 
the building blocks (molecules and granulates) are sim- 
ilar, there are other important differences between the 
two systems, e.g., thermal fluctuations are important for 
molecules while they are negligible in the granular case. 
Also, the interparticle forces have very different charac- 
ter: the molecules experience attractive forces (dipolc- 
dipole, van der Waals), while the granular particles in- 
vestigated here interact only via hard core repulsion and 
friction. 

As seen in Fig. [5j the average alignment angle av de- 
pends only on the properties of the particles. The shear 
rate independent av has been determined for all eight 
samples and it is presented as a function of the aspect ra- 
tio L/d in Fig. [6] The plot combines the bulk and surface 
data from x-ray tomography, the surface data from opti- 
cal experiments, and the numerical results in bulk plane 
Couette flow. We find a systematic decrease of 6 av with 
increasing length to width ratio for a given shape family 
(rice or cylinders). Results of the numerical simulations 
nicely agree with the experimental data. At comparable 
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FIG. 7: (color online). Evolution of the particle orientations 
as a function of strain 7 in the stationary state for glass 
cylinders with L/d — 3.5. Data obtained by optical particle 
detection. The black line corresponds to dO/dry = —1 (in 
radians). For better visualization of the trajectories we used 
4 different colors. 



the orientational distributions. This range is seen as a 
dense band in Fig. [7] between the angles of —20° and 
35° for the case of glass cylinders with L/d = 3.5. In 
this range the particles rotate forward or backward with 
essentially equal probability. Outside this range the par- 
ticles are rarely seen and they rotate relatively quickly 
forward as it is indicated by the steep trajectories in 
Fig. [7J For these measurements data have been taken 
from the center of the shear zone, where the shear rate 7 
is nearly constant. 
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FIG. 6: (color online). The average orientational angle av as 
a function of the length to diameter ratio L/d of the particles 
obtained by experiments (»,A,A) and numerical simulations 
(x). 

aspect ratios, the alignment angle appears to be system- 
atically smaller for cylinders than for ellipsoids. 

Focusing on the dynamics of individual grains during 
stationary shear, the orientation of the majority of the 
particles stays within a range of O av — a < 6 < av + a, 
where a is the standard deviation of the Gaussian fit of 



FIG. 8: Distribution of the rotation velocity to of the par- 
ticles as a function of their orientation 9 in the stationary 
shear aligned state for glass cylinders with L/d — 3.5. Data 
obtained by optical particle detection. In the grayscale plot 
black corresponds to the highest probability density. 

The distribution of the particle rotation velocity u) = 
&6/dt has been determined as a function of and the 
resulting plot is shown in Fig. [8] Here cj is normalized by 
the average shear rate at the location of the particle. As 
it is seen ui/j is better defined for the above mentioned 
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range (9 av — a < 9 < 9 av + a) near av , while the 
distribution is wider (i.e. uj is less well defined) near 
= ©av + 7r/2. At 6 = 8 av + 7r/2 we find an average 
rotation speed cj/7 w —1 rad/unit strain. Thus, at this 
angle the rotation speed of the particles is maximal, it 
reaches the value corresponding to the local vorticity of 
the shear flow. 

We calculated the averages of the rotation velocity dis- 
tributions normalized by the shear rate, w av /7, for var- 
ious materials, and plotted them against 9. Data from 
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FIG. 9: (color online). Normalized average rotation velocity 
of the particles i^av/7 as a function of their orientation 6 for 

(a) rice and glass cylinders obtained by optical measurements 

(b) pegs obtained by x-ray tomography, (c) numerical simula- 
tions and (d) hard ellipsoids in a sheared liquid as calculated 
by Jefferv [25]. 



optical measurements are presented in Fig. [9ja) for rice 
grains with L/d = 2.0, 3.4, 4.5, and glass cylinders with 
3.5, while Fig. [9^b) shows the results obtained by x-ray 
tomography for pegs with L/d = 2.0, 3.3 and 5.0. Fig- 
ure pnc) presents results obtained by numerical simula- 
tion (DEM). We have also included the result of the cal- 
culation by Jeffery [29] , which describes the rotation ve- 
locity of a hard prolate ellipsoid immersed in a sheared 
viscous liquid. All four panels of Fig. [9] show similar ten- 
dencies, viz. the particles rotate faster when their ori- 
entation is far from the preferred one and rotate slower 
when they are closer to it. The ratio of this fast and slow 
rotation speeds appears to be smaller for shorter grains. 
In other words, longer grains show a stronger angular de- 
pendence of the rotation rates. Differences between our 
results for the granular particles and the Jeffery solutions 
arise from the grain-grain interactions. 



B. Transients into the aligned steady state 

So far, we have discussed the properties of the sta- 
tionary state. In the following, we focus on the dynam- 
ics observed in the transient during which the system 
reaches the above described stationary state. Two types 
of experiments have been carried out: (a) starting from 
an initially random configuration, and (b) reversing the 
shear direction, i.e., starting from an oppositely shear 
aligned state. The evolution of the distribution of parti- 
cle angles has been investigated for glass cylinders with 
L/d = 3.5 taking (a) 1350 and (b) 1700 independent 
measurements, respectively. These measurements were 
performed using surface optical imaging, and the result- 
ing probability distributions are plotted on Fig. [lO] as a 
function of the strain (similarly to Fig. J7|. 
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FIG. 10: Evolution of the orientational distribution of the 
particles (a) starting from an initially random orientation (b) 
from an aligned state after shear reversal. Both measurements 
were performed with glass cylinders with L/d = 3.5, data ob- 
tained by optical particle detection. In the grayscale plot, 
black corresponds to the highest probability density. The 
white arrow indicates the evolution of the center of the dis- 
tribution for case (a). For case (b) the initially ordered state 
appears to split into two independently relaxing distributions, 
marked by white arrows. 

For the initially random system (a) we find a continu- 
ously increasing order and a decreasing alignment angle 
[see white arrow in Fig. |10[a)] both converging towards 
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their stationary values. The evolution of the order pa- 
rameter S and that of the average angle 9 av are shown 
as functions of 7 in Figs. [Tl|a)-(b). As it is seen, or- 
dering develops faster than the convergence of the shear 
alignment angle. This can be quantified by the character- 
istic strain parameters of the approximately exponential 
relaxation, which are 75 = 0.73 and 79 = 2.57, respec- 
tively. 

For initial condition (b) we find that the average angle 
rotates in the opposite direction, from — O av to av - he. 
the majority of the grains turns backwards and only a 
smaller fraction rotates forward as indicated by the two 
white arrows in Fig. 
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During this procedure the 
order parameter initially decreases then increases back 
to its stationary value as seen in Fig. 11 'a). Interestingly 
the strain which is necessary to reach the stationary an- 
gle is about the same for the two initial conditions, but 
the strain needed to reach the stationary value of the or- 
der parameter is significantly different. To destroy and 
then recover the order takes about three times as much 
strain (75 = 2.19) compared to alignment from an ini- 
tially random state (75 = 0.73). 
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FIG. 11: (color online). Evolution of the (a) order parame- 
ter and the (b) average alignment angle during both type of 

(c) Evolution of 
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transients. Same measurement as in Fig. 
the applied torque needed to maintain stationary rotation of 
the bottom plate starting from an initially random orienta- 
tion. All measurements are for glass cylinders with L/d — 3.5. 
Continous lines are fits assuming exponential relaxation, see 
text for the numerical values of the exponents. 



It is very interesting to compare the normalized aver- 
age particle rotation speed w av /7 as a function of 9 for 
these two different transients and the steady state, see 
Figs. |12[a)-(b). In the stationary state, grains in all ori- 



entations 9 rotate on average with negative u> (cf. Fig. 12 
black circles). During the transients, we find a positive 
mean rotation velocity for grains between —15° and 15°. 
This means, that for both cases (a) and (b) these par- 
ticles rotate against the direction dictated by the shear 
flow, while the rest of the particles rotate forward. Here, 
for the transients, the difference between cases (a) and 
(b) is that for the "oppositely shear aligned" transient the 
number of backward rotating grains is much larger due 
to the nature of the initial configuration. This is why the 
average reorientation of the grains 9 av ("f) differs in the 
two cases (see white arrows indicating the grain rotation 
direction in Fig. 10 and the curves in[Tl|b)). 
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FIG. 12: (color online). Normalized average rotation veloc- 
ity of the particles w a v/7 as a function of their orientation 9 
during the two types of transients (a) starting from a random 
orientation (b) starting from a reversed flow aligned state. 
Measurements done for glass cylinders with L/d = 3.5, data 
obtained by optical particle detection. 

Additional measurements were performed for the ini- 
tially random configuration (a). In these measurements 
we aimed to monitor the evolution of the effective friction 
of the sample during the transient. For this purpose we 
measured the torque M needed for a stationary rotation 
of the bottom plate, which was rescaled by the asymp- 
totic value Moo as M = M/M^. A reference dataset was 
also recorded for nearly spherical mustard seeds, which 
is used to correct for effects not directly related to the 
reoriantation of the sample. The rescaled and corrected 
curve is shown in Fig. 
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c). The first thing to note 
is that the ratio of the initial and final values of M is 
around 1.5, suggesting that the effective friction of the 
unoriented sample is about 50% higher than that of the 
shear aligned one. The second observation is, that the 
evolution of the effective friction of the material can be 
described by an exponential decay of the excess torque as 
M — Mmustard « (M max - l) e -^ M . This exponential fit 
results 7m = 1-46 which is in between the characteristic 
values obtained for 73 and 79. Altough the experimental 
curve is fairly noisy, one can pursue a more detailed anal- 
ysis by fitting the data with the sum of two exponentials 
with decay rates of 7s and 79. The resulting amplitudes 
suggest that the evolution of S has a stronger influence 
on the evolution of M compared to the contribution of 



V. CONCLUSIONS 

The orientation and rotation of elongated granular par- 
ticles has been tracked in shear flow. In the asymptotic 
state, the time and ensemble averaged orientation of the 
particles forms a small angle O av with the streamlines. 
This shear alignment angle is independent of shear rate 
and decreases with increasing particle aspect ratio. The 
characteristics of the asymptotic alignment determined 
from optical experiments in this study are in good agree- 
ment with experiments obtained from x-ray computed 
tomography |30| of cylindrical pegs. 

The average rotation speed of individual particles as 
function of the orientation uj av (9) shows that particles 
rotate faster if their orientation is perpendicular to the 
preferred average orientation and their rotation is slower 
when their orientation is near O av . The ratio of this 
fast/slow rotation increases with grain aspect ratio, sim- 
ilarly to the case of a rotating hard ellipsoid in a sheared 
liquid as calculated by Jeffery [29 . In the stationary 
shear aligned state the net rotation rate of particles is 
always negative i.e. the sense of rotation is dictated by 
the shear flow. Transients into the asymptotic alignment 



have been studied starting from an (a) initially random 
configuration and (b) an oppositely shear aligned state. 
The character of the transients depends upon the initial 
order and alignment. The transient is longer for case (b) 
where initial order is partially destroyed and then recov- 
ered again. A common feature of the two transients is 
that the particles with orientation just beyond av rotate 
backwards, and reach the preferred angle on a shorter 
way. 

The evolution of the mechanical resistance of the sam- 
ple measured during the transient indicates a 30% reduc- 
tion of the materials effective friction during the course 
of the shear induced ordering process. 
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